Entanglement and Classical Polarization States 
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We identify classical light fields as physical examples of non-quantum entanglement. A natural 
measure of degree of polarization emerges from this identification, and we discuss its systematic 
application to any optical field, whether beam-like or not. 
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From Christian Huygens' explanation of the fascinating 
birefringent property of the crystals called Iceland spar 
(i.e., calcite), through the much later work of Sir George 
Stokes [l| , the formulation of polarization theory has con- 
tinuously evolved and is now well established in terms of 
field correlation functions @j. The concept of degree of 
polarization is based on the coherence matrix or polariza- 
tion matrix constructed from these functions. But, even 
after centuries of attention to this basic property of opti- 
cal fields, fascinating new issues concerning polarization 
have emerged in the past two decades. 

The familiar measures of polarization come from the 
treatment of light as a beam. This implies a given direc- 
tion of propagation, and thus a specific transverse plane. 
But development of highly non-paraxial fields, use of very 
narrow-aperture imaging systems, recognition of associ- 
ated propagation questions, and probing of fully three- 
dimensional fields as in hohlraums, all point to necessary 
modifications of polarimetry beyond the traditional pic- 
ture [3|— 15|] - 

Questions of definition and of principle are also in need 
of answers. For an electromagnetic field without a clear 
transverse plane or even cylindrical symmetry the rel- 
evance of measures such as the conventional degree of 
polarization must be reconsideredjff , and alternative ap- 
proaches have been proposed |7H12| . The same is true of 
experimentally important measures, such as the Mueller 
matrices [l3|, [l4[ • In addition, pioneering measurements 
[l5l - [l8j have highlighted two-photon and multi-photonic 
views of polarization (l9l [2pj , and polarization measures 
for nonlinear classical optical waves have been considered 



2 1| | . Here we describe a basis for comprehensive consid- 
eration of interrelated questions by addressing the issue 
of dimensionality in a unified way. 

Of course, dimensionality is trivially engaged in con- 
verting the field from planar-transverse to non-planar, 
which only requires the natural extension 



E = xE r 



yE y 



xE x + yE y + zE z 



(1) 



in order to take into account a third component of the 
field. However, an entirely comprehensive treatment of 
polarization questions should begin by first noting that 
two independent vector spaces are employed in each re- 
alization of E, and second that the two spaces are en- 



tangled [22j. Entanglement as a technical term means 
just that E, in (JXJ) , is a tensor product of "lab space" 
unit vectors such as x and y, and functions E x and E y 
that are vectors in a statistical "function space" of con- 
tinuous normed functions (typically taken to depend on 
space-time or space-frequency). 

The light field uses at least one vector from each of 
these distinct spaces, and in the general case it is not 
possible to convert ([T]) to the completely polarized form 
E = uF, in which the two spaces appear only in a fac- 
tored direct product. Such a (typically unlikely) form 
is called by the equivalent terms "separable" or "fac- 
torable" or "non-entangled" . That is, determining degree 
of polarization is the same as determining degree of fac- 
torization (separability) of the two spaces, i.e., absence 
of entanglement. This fact can be exploited. 

Both inner and outer products play a role in polariza- 
tion theory. For convenience we first consider a beam 
and write the fundamental polarization outer product, 



|E)(E| 



x)|£L 



;x|(^| + ( y |(^| , (2) 



as the field intensity times a normalized hermitian outer 
product W, i.e., I W = |E)(E|. Here the angle brackets 
explicitly make the point mentioned above that the unit 
vectors and the field components are members of different 
vector spaces. This is a quantum-like notation that will 
be helpful, but no quantum properties will be introduced. 

With the intensity I = (E x \E x ) + (Ey\E y ) factored out, 
we can write 



V7 cos#|x)|ea;) + sin0|y)|e 



IE) 



and this allows W to be written: 



(3) 



W = fcos0|x)|e x ) +sin0|y)|e„)V<x|(e a! |cos0 



{y\{e y sine/ , 



(4) 



where assignment of the relative amplitudes via sine and 
cosine factors allows the components |e x ), \e y ) to be unit- 
normalized: (ej|ej) = 1. We take account of the generally 
non-zero correlation between the field's components by 
introducing the magnitude and phase of the cross corre- 
lation as 



= a = \a e 



(5) 
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Although rarely mentioned, each of the two separate 
vector spaces has its own polarization matrix. These are 
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reduced-state tensors, i.e., traced over one space inde- 
pendent of the other. The "normal" polarization matrix 
is obtained by tracing over the function space. We can 
denote it as Wi a b = 7Yf cn (W), and calculate it by a diag- 
onal sum over any complete set of orthonormal vectors in 
the function space (but only in the function space), a set 
that we can label {|0 m )}. For short we will temporarily 
use p and q to stand for x or y, and then obtain 

|E)(E|0 ml 

= V Piq \p){q\{Eq\E p ), {p,q = x or y), (6) 

which we recognize as a 2 x 2 tensor in the lab space in 
the basis defined by |x) and |y). When the intensity is 
factored out we have a familiar matrix expression with 
the required unit trace: 



W lab = Tr fcn (W) 



cos a 
a* sin 9 cos ( 



a cos o sm f 
sin 2 9 



(7) 



The less familiar other polarization matrix is obtained 
by tracing over the lab space. This is trivially done via 
the projections |x)(x| and |y)(y|, but the result will not 
be in standard form because of the non-correlation of E x 
and E y . We can overcome this by rewriting |E) as a sum 
of a pair of statistically orthogonal components. If we 
choose \e x ) as one of them, we will denote its partner by 
\e x ), with (e x \e x ) = 0. Then \e y ) becomes a combination 
of both components, a\e x ) + f3\e x ), so 

|E)/V7= (cos#|x) +asin0|y))|e a: ) + /?sin%)|e x ), (8) 

where a is defined in (JSJ) , and \a\ 2 + \/3\ 2 = 1. Then 
the reduced polarization tensor for the function space 
Wfcn = T>iab(W) has the form: 



W fcn 



COS 



|a| 2 sin 2 6> a*/3sin 2 



a/3* s\n z 9 



| 2 sin 2 



(9) 



These results serve as background in addressing polar- 
ization of non-paraxial or non-beam light fields, i.e., the 
three dimensional fields expressed as given in (fTJ), and 
non-trivially entangled. We recall the observation that 
(TT]) has the character of an entangled description of two 
parties jointly, where the "parties" here are the two dis- 
tinct aspects (or degrees of freedom) of the light field, 
namely the lab space direction of the optical field and 
the statistical function space characterizing the strength 
of the field. This suggests making a Schmidt-type analy- 
sis [231 124| , as frequently used in discussions of quantum 
entanglement data in few-mode and multi-mode photonic 
contexts Hi-Hi]. 

The Schmidt theorem applies to any kind of two-party 
vector, whether quantum or not, and we begin by calcu- 
lating the same two reduced matrices Wi a b and Wf cn de- 
fined above, except that now they are three-dimensional. 
Since they arise from a common hermitean W, which 
is the 3x3 analog of (j3|, they share the same three 
(real) eigenvalues, k 2 ,k 2 ,k 2 . Their eigenvectors [3(| are 



orthonormal and occur in pairs. The Schmidt-form re- 
sult enables the optical field to be written immediately in 
a perfectly organized form in terms of these eigenvalues 
and eigenvectors: 



IE 



Kl|U!)|/i) + « 2 |u 2 )|/ 2 ) + K 3 |u 3 >|/ 3 ). (10) 



Here perfect organization means that orthogonality con- 
ditions apply in both vector spaces at the same time: 
(fi\fj) = (uj|u.j) = Sij, and since intensity has been fac- 
tored out, the three ks are normalized on the surface of 
a unit sphere: ftf + k 2 , + K 3 = !• 

This Schmidt decomposition of the field is unique up 
to a rotation at most, and allows interesting questions 
to be answered by inspection. For example, it is obvious 
that (JTUJ) can take the completely polarized (fully fac- 
tored and so non-entangled) form E = uF = \fl u / only 
when two of the ks have the value zero. Because the three 
|/) vectors are mutually statistical orthogonal, no other 
k values can produce a completely polarized result. For 
such a field, no matter which lab direction v is used for 
a projective measurement of E, the statistical-functional 
features of that component will always be exactly those 
of the function /. And similarly, no matter what projec- 
tion in function space is employed, the projected field's 
direction will always be u in lab space. 

The connection to standard polarization measures is 
not complicated. Obviously the conventional definition 
2] of degree of polarization P, designed for a beam-type 
field, is sensible only when there is zero projection along 
one of the three it's, say along w 3 , in which case k 3 = 0. 
Then P has several compact expressions @, [H, |3l| : 



P 2 = l-4Det(VVU) 



2 1 2 



| 2 = 1-2(1-1/X), (11) 



where K will be introduced below. As is well known 
0, the beam-based definition resists being generalized 
to the three dimensional case when none of the KjS is 
zero in (|10[) . In addition, with k 3 = 0, the field cannot 
point even slightly into the direction u 3 . Thus to call 
it "unpolarized" is not fully sensible even if K\ = k 2 , 
making P zero. 

The Schmidt decomposition automatically provides a 
very useful "weight" parameter K [32[ , which counts the 
non-integer effective number of dimensions needed by the 
optical field. The expression for K is 



K = 1/[kI + i4 + k 



A] 
3J> 



(12) 



which gives greater weight to the vector directions with 
the larger absolute K values. It's easy to see that K lies 
between 1 and 3 and incorporates the beam-type field 
case automatically when any one of the three k's is zero. 

K = 1 occurs if k 2 = k 3 = 0, which signals a one-term 
E, i.e., completely polarized light, obtainable only if the 
original field components could have been rotated into a 
fully factored form: E = uF. This is of course generally 
not possible. The farthest departure occurs for the value 
K = 3, when k 2 — 1/3 for all three components, and 
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FIG. 1: The entanglement measure K varies from 1 to 3 over 
the unit polarization sphere. The purple zone-centers touch 
the surfaces of the cube (K=l, completely polarized), and 
the red centers are completely unpolarized (K=3). Planar- 
polarized fields (K = 2 or P—Q) are located at the corners 
of the triangular regions where one Kj = and the others 
equal 1/2. The mesh lines locate partially polarized values 
K = 3/2,2,5/2. 



all fj are equally intense. This is maximal entanglement 
and also what is sensibly called a completely unpolarized 
field. Intermediate values of K represent intermediate 
degrees of entanglement (partially polarized light). All of 
these conditions are associated with k values that identify 
points on the unit sphere in Fig. [T] 

In summary we have reformulated polarization the- 
ory as entanglement analysis. The Schmidt theorem ap- 
proach automatically provides an optimum expression for 
any light field by identifying its orthogonal directions in 
lab space, u^, and its associated amplitudes ft in statisti- 
cal function space. This is possible because every optical 
field is a quantity existing simultaneously in those two 
independent vector spaces. The interpretation of degree 
of polarization naturally corresponds to the separability 
between the two spaces for both planar and non-planar 
cases. In this new perspective, polarization is a charac- 
terization of the correlation between the vector nature 
and the statistical nature of the light field. 
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